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Abstract. The Perron-Frobenius theorem is apphed to identify the superfluid 



^ ' transition of a two-component Fermi gas with a zero-range s-wave interaction. 



According to the quantum cluster expansion method of Lee and Yang, the grand 
partition function is expressed by the Lee- Yang contracted 0- graphs. A singularity of 
an infinite series of ladder-type Lee- Yang contracted 0-graphs is analyzed. We point 
out that the singularity is governed by the Perron-Frobenius eigenvalue of a certain 
primitive matrix which is defined in terms of the two-body cluster functions and the 
Fermi distribution functions. As a consequence, it is found that there exists a unique 
fugacity at the phase transition point, which implies that there is no fragmentation 
of Bose-Einstein condensates of dimers and Cooper pairs at the ladder-approximation 
level of Lee- Yang contracted 0-graphs. An application to a Bose-Einstein condensate 
of strongly bounded dimers is also made. 
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1. INTRODUCTION 

When physically relevant operators are primitive in some representation, the Perron- 
Frobenius theorem may play a role in the study of physics. Examples include theories 
of stochastic processes P and studies of ground states of quantum many-body systems 
[2],[3]. In this study, we provide a specific example of primitive matrices, which belongs 
to a different category. It concerns a two-body cluster function in quantum-statistical 
mechanics, which is important to dilute Bose and Fermi systems. 

To be specific, we consider two-component atomic Fermi gases with a zero-range 
interaction characterized by the s-wave scattering length a. This system exhibits a 
crossover between a BCS-like state of Cooper pairs and a Bose-Einstein condensate 
(BEC) of diatomic molecules [1]. The two-body cluster function U^"^^ mentioned 
above appears in the quantum cluster expansion of the grand partition function. The 
quantum cluster expansion is important to investigate thermodynamic properties of 
the system, and it has seen a resurgence of interest recently [51-fTO]. The cluster 
integrals h^'^^ and 6*-^^ for the unitary Fermi gas have been experimentally measured 
using a mixture of ^Li in two internal states [HI [12]. Here the cluster expansion is 

Pp\^ = 2+ '^z^ + b'^^h^ + b^^h"^ H , where p is the pressure and A := {2Tr/mkBTy^^ 

is the thermal de Broglie length. The second cluster integral b^'^^ = can be calculated 
from the trace of the two-body cluster function ?7*-^^ . In this regard, the two-body cluster 
function U^'^^ has more information than the second cluster integral b^'^\ Under the s- 
wave approximation, we show that t/*-^^ is a non-negative matrix with the Bethe-Peierls 
boundary condition [13] or the Fermi pseudopotential [Ti] . 

According to the quantum cluster expansion method of Lee and Yang [T51IT6] . 
the grand partition function is calculated by the cluster functions U^^\ U^'^\ etc., 
and expressed by the Lee- Yang contracted 0-graphs. In this study, we consider an 
infinite series of the ladder-type Lee- Yang contracted 0-graphs, which gives the non- 
interacting dimer BEC transition temperature in BEC limit [I71[T8], and the BCS- 
transition temperature in BCS limit [I9]. To investigate the phase transition of the 
system, we analyze a singularity of an infinite series of ladder-type Lee- Yang contracted 
0-graphs. We point out that the singularity is governed by the Perron-Frobenius 
eigenvalue of a certain primitive matrix which is defined in terms of the two-body cluster 
functions U^'^'^ and the Fermi distribution functions np. As a consequence, there exists 
one and only one fugacity at a phase transition point and physically it implies that 
there is no fragmentation of Bose-Einstein condensates of dimers and Cooper pairs at 
the ladder-approximation level of the Lee- Yang contracted 0-graphs. 

This paper is organized as follows. In section 2, the phase transition is defined 
as the disappearance of holomorphy in thermodynamic functions. In section 3, the 
Lee- Yang cluster expansion method is reviewed for the case of two-component Fermi 
systems. In section 4, the matrix elements of the two-body cluster functions in 
momentum representation are calculated and are found to be non-negative. In section 
5, a singularity of an infinite series of the ladder-type Lee- Yang contracted graphs are 
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analyzed by using the Perron- Frobenius theorem. An apphcation to the BEC of strongly 
bounded dimers is also made. In section 6, we summarize the main results of this paper. 
The proofs of several formulas and inequalities are relegated to appendices to avoid 
digressing from the main subject. 



2. Grand partition function and phase transition 

To be specific, we consider a system of two-component Fermi particles with the same 
mass m in a finite volume = at a temperature T. We also introduce the fugacity 
z = e^^, where /3 = l/Zc^T is the inverse temperature and n is the chemical potential per 
particle which is assumed to be independent of spin variables a =t", Our A^-particle 
Hamiltonian and the grand partition function are given by 

^^"^ = -^Ev^ + E-(ir^-r.i), (1) 



Sv^:=l + ^.^Tr^(.,[e-'^^'^'], (2) 

N=l 

where Tr (jv)[- ■ ■] denotes the trace over the antisymmetrized A^-particle Hilbert space 

7/^\ Here we set /?. = 1. The equilibrium pressure p and the particle-number density 
p of the system are given by 

P = Jim ^log^^' (3) 

and 

p = lim ^z-^hgEy. (4) 

v-^oo V dz 

A phase transition is mathematically characterized by singularities in thermody- 
namic functions. The singularities here imply the disappearance of holomorphy. A 
complex function defined on C is holomorphic on an open set f2, if its Taylor expansion 
around every point in VL has a nonzero radius of convergence. If Vl is not an open set 
we interpret that holomorphy holds in an appropriate open set containing VL. In lit- 
erature, the term 'analytic' is often used instead of holomorphic; however, we use the 
latter in this paper, because the former is a global concept. In the classical limit, the 
cluster expansion of p has a positive radius of convergence [20], so if the fugacity z is 
sufficiently small, p is holomorphic in z. That is, there is no phase transition if 2; ^ 1. 
If a singularity in a thermodynamic function such as p{z) exists at a point Zq on the 
positive real axis and the holomorphy disappears there, we conclude that there is a 
phase transition. Lee and Yang analyzed zeros of 5 for classical lattice gas systems (or 
equivalently classical Ising systems [21]. In contrast, for the superfluid phase transition, 
we analyze zeros of 

Now we assume that a binary interaction v is characterized by the s-wave scattering 
length a and that all other relevant length scales, such as the thermal de Broglie length 
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A := ^j2Tiffi/mT and the mean interparticle distance, are much larger than the range 
of interaction. Then, by dimensional analysis, we have 

lim ^ log = (z. A/ a) . (5) 

The superfluid transition manifests itself through the singularity of f{z,X/a) on the 
positive real axis of z. 

3. Lee- Yang cluster expansion method for a two-component Fermi system 

In the Lee- Yang cluster expansion method, the grand partition function for a system 
of Fermi particles can be expressed in terms of the contracted 0-graphs [16]. Each 
contracted 0-graphs is computed from the Fermi distribution function np{k.) := {1 + 
^-ig^k2/(2m,)^-i A^-particle cluster functions U^'^^ for the same system obeying 

Boltzmann statistics. 

In this section, we first define the A^-particle cluster functions U^^\ Then, we show 
how to express the grand partition function in terms of the contracted 0-graphs. 

3.1. Definition of Cluster functions 

To define the A^-particle cluster function, we first introduce 

(1, . . . , iV|iy(^) |1, . . . , iV) := iV! 5^ (!', ...,N') ^* (1, . . . , iV) e-^^% (6) 

where we use 1,2,... and 1', 2', . • • to represent the momentum and spin coordinates of 
individual, 1 := (ki,(Ti), etc. and 1' := (k[,a[), etc., and ipi and Ei are the normalized 
eigenf unctions and eigenvalues of H^^^ , respectively. The summation in Eq. ([6]) extends 
over all eigenf unctions in an unsymmetrized iV-particle Hilbert space 7i^^\ The 
matrix elements of U^'^'^ in momentum representation are 

(l'|f/W|l) := (l'|e-^^'"|l), (7) 
(l',2'|f/(2)|l,2) := (l',2'|e-''^'"|l,2) - (l'|t/«|l) (2'|f/«|2), (8) 
(l',2',3'|t/(3)|l,2,3) := (l',2',3'|e-'^^^"|l,2,3) - (l'|f/W|l)(2', 3'|t/(2)|2, 3) 

- (2'|f/(i)|2)(3',l'|f/(2)|3,l) - (3'|f/(^)|3)(l',2'|t/(2)|l,2) (9) 

- (l'|t/«|l)(2'|t/W|2)(3'|t/W|3), etc.. 

In the computation of the grand partition function for Fermi particles, only the anti- 
symmetric combination T^'' appears. We define the matrix elements of T^'' by 

(1', . . . , r |Tf ) 11,...,/):=^ {-lf{P{l'), • • • , Pin |f/^^^ I 1, . . . , / ), (10) 

P&Sn 

where P denotes permutation and (—1)^ = 1 or —1 for even or odd permutations P. 
We discuss a couple of examples. 



Example 1 — 
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The one-particle cluster function is 

(kV|f/(i)|ka) = (5.v(k'|t/(i)|k) = 4,k'5.,.'e^^^'/(2'"). (11) 
We consider the geometric series 

oo -. 
^(^_^g-/3kV(2Hj"^ (^l^^^-/3kV(2r„)j- =l_^p(k). (12) 

n=0 

The effect of Fermi statistics emerges through this geometric series. 

Example 2 — 

The momentum part of the two-body cluster function is given by 

(k;(T;,k'2a^|t/(2)|kiai,k2a2) = 5.,,, 5,,^ (k^, k'^jf/^^) |ki, ka). (13) 

In this paper, we discuss the case of two-component Fermi particles with an s- 
wave approximation. Then, a two-body cluster function is a symmetry matrix 
(k'i,k'2|f/(2)|ki,k2) = (k'2,k;|f/(2)|ki,k2), and 

(l',2'|Ti')|l,2) = (k;,k'2|f/(2)|ki,k2) ■ ((5.,.i<5.,.^ - . (14) 

The two-particle Hamiltonian is 

Hi^) = + + vir) = + H^'"''^ (15) 

i/M) = lv2 + ,;(r). (16) 
m 

The two-body cluster function for a finite volume V is 

(k;k'2|[/(2)|kik2) = ^5K,K'e-^^^/(^™) • (k>(-i)|k), (17) 

(k>(^^')|k) := (k'le-^^'^'^lk) - 5k,k'e-^'''/'", (18) 

where K := ki + k2 and k := (ki — k2)/2, and the Kronecker delta reflects the 
conservation of momentum. The function (k'|u''''^^)|k) describes an effect of interaction 
and can be calculated from the eigenf unctions and eigenvalues of H^^'^^^ [T3] . 

3.2. The grand partition function represented by contracted 0-graphs 

In terms of the cluster functions, the grand partition function S is written as [T5|[T6] 

logSy = -2^Li5(-^)+P, (19) 

where Li;(a;) := 'Yl^=i^^ /^^ is the polylogarithm and the first term on the right-hand 
side gives the logarithms of the grand partition function of a non-interacting Fermi gas. 
Here V is the sum of all different contracted 0-graphs introduced by Lee and Yang [16] 
and is illustrated in figure [TJ The number under each graph in figure [1] shows the 
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Figure 1. (a) Sum over all different contracted 0-graphs (dotted curves). The 
symmetry numbers are shown under the graphs, (b) Expression of (k'(T'|7/o|kcr) as the 
sum over different primary 1-graphs (solid lines). 



symmetry number of the corresponding contracted 0-graph. The algebraic expression 
of the sum of the contracted graphs is 

^ = y ^ ^ ?7o(ki)?7o(k2)(kicri,k2cr2|TA^|kicri,k2cr2) 

ki,k2 a-i,(T2 

+ "g" 5Z X] ?7o(ki)?7o(k2)?7o(k3)(ki(Ti,k2a2,k30-3|TSJ^|kiCTi,k2cr2,k30-3) 

ki,k2,k3 0-1,0-2,0-3 

+ "2" ■ ■ ■'7o(k4)(kicri,k2cr2|T^^|kicri,k3a3)(k3(T3,k4cr4|T^^|k2cr2,k4cr4) 

ki,...,k4 0-1,..., 0-4 

+ y ^ ^?7o(ki) ■ ■ ■77o(k4)(kiCri,k2a2|T^^|k3Cr3,k4Cr4)(k3Cr3,k4Cr4|T^^|kiCri,k2Cr2) 
ki,...,k4 0-1,. ..,0-4 

+ ■■■, (20) 

where each term in the sum corresponds to the contracted 0-graph at the corresponding 
order in figure [H Here, 

r/o(k) := 1 -nF(k) = (I + z^-^''"'^^^^)-^ (21) 

describes the effect of the Fermi statistics. 



4. Two-body cluster function in the s-wave approximation 



To calculate an analytical representation of (k'^, k^f/'^^^lki, k2) within an s-wave 
approximation, we can use the Bethe-Peierls boundary condition [13] 



1 d 

rip dr 



r-s.+0 



1 

a 



or the Fermi pseudopotential 



, , 47ra ^3 , , d 

Vps{r) := (r)— r, 

m or 



(22) 



(23) 
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Unitary , 




BCS 









-2-10 1 2 



w 

Figure 2. The logarithms of three terms in the function s{x^w). "BCS", "Unitary", 
and "BEC" show — log2, — log(|w|) + 0.3624, and + log(ierfc(— w)), respectively. 
Here 0.3624 ~ yisyi{xF [x) j ypH) {x ~ 1.502). The dashed curves are - log(2V7r |w|) 
and up' . 



where a is the s-wave scattering length. The Bethe-Peierls boundary condition or the 
Fermi pseudopotential supports continuous scattering states 

V^sc('^) = {27r^[l + (/csco)^]} 2 ^(sin(/cscr) - /Cscacos(fcscr)), (24) 

with the energy E^c = k'^^/^, and if > one bound state 

^b(r) = (27ra)-^-e"^/^ (25) 
r 

with the binding energy E-^ = — l/(ma^). Using the set of energy eigenvalues and 
eigenstates ( l24l) and ( l25l) . the two-body cluster function can be computed as 

(KXluS^l^iM) = ^5K,K'e-^^^/^ • (k'|4f |k), (26) 
lk\ - s{x',w)-s{x,w) 

^ w e"^' - ^xF{x) -w e^-'erfc {-w)] . (28) 



where 



2;2 _|_ , 



Here we introduce the dimensionless variables x := ^J'^WJm = Xk/y/Ojr and w := 
\fP) I i^^fma) = A/(v^27ra), and the Dawson's integral 

F(x) =e-"' r dte'\ (29) 







and the complementary error function 



erfc(x) = 4= / dt^~*^- (30) 



The proof of fl2Bl) - fl251) is given in Appendix A The two-body cluster function in the 
s-wave approximation has been first obtained in Ref. [17] for > and in Ref. [22] 
for < 0. However the representation in f l26p -f l28p has two advantages over that 
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CO 



CO 




Figure 3. Dependence of -g^s{x,'w) on x := ^ fik'^ /m = Afc/\/27r for w := 
\/^/{^/Tna) = \/{\/2TTa) = —4, —3.5, —3, . . . , 1, 1.5. The dashed curve shows the case 
of the unitarity hmit {w ~ 0). We see that ^s(a;, w) > for arbitrary — oo < w and 

< X < GO. 



in Ref. [T7tl22] : (i) it holds for either sign of a and (ii) it can easily be seen that 
the two-body cluster function is continuously connected at the unitary limit = 0. 
The function s{x,w) in ( 128|) is composed of the three terms. In the BCS region 



(A/a = y/2n (3 / (y/ma) < —1), the first term (i.e., e '^^) is dominant. In the unitarity 
region (A/ \a\ < 1), the second term (i.e., ~^^F (x)) is dominant. In the BEC region 
(A/a > 1), the last term (i.e., — e'^^erfc (— w)) is dominant. The logarithms of these 
three terms are shown in figure [2l 

The derivative ^s{x, w) for various values of w is shown in figure [3] which indicates 
that -^_s{x,w) > except at x = 0, and hence s{x,w) is a monotonically increasing 



function of x. In fact, this can be proved [23] and is shown in [Appendix B[ Moreover, 



we have ^s(0,w) = 0. Thus, the following theorem is established: 



Theorem 1.- 



For fixed u; > — oo and any k, k' G (k'lujj'^lk) > 0, with the equality holding 

if and only if k = k' = 0. 

Theorem 1 plays a pivotal role in the main result of this paper in Sec. [5l 
It should be noted that Theorem 1 is established by using the exact energy 
eigenvalues and eigenf unctions. Therefore if some crude approximations are employed, 
the above theorem may not be derived correctly. For example, at the first-order 
perturbation in a, (k'|Mgp^|k) is not a non-negative matrix [23]. To see this, we first 
calculate 

{KK\5\r, - r2)|:r|ki,k2) = ^2j/^^j/^^ e^(^-^'>^e"^'-5^(ri - v,)^^ [r^^^) 
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= :^5K,K'A(k). (31) 

Here A(k) is the Tan's A-function, which is defined as A(k) = 1 for A; < oo, together with 
j (i^k/c~^A(k) = [25]. Clearly, A(k) becomes negative for sufficiently large k. In fact, 
the Tan's A-function is written as A(k) := 1 — k^^5{k^^) [26]. We denote (k'|M^gf^|k) 
by (k'|Mgp''|k) = (k'lM^gf-* |k) + O(a^). Using the general formula discussed in Ref. [15] 
and equation fl3T|) . we obtain 

(k;,k;if/S|ki,k2) 

= - f rfr e-(^--)«'+'=")/(2™)^]^/ ^ KM^i - r2)|ki, ks) 
Jo 

^''-5K,K'e-^^^/(^-)-(k>£f|k), (32) 



V 



where 



(k>Ss?|k) = ^ ^ A(k). (33) 



27r^m k''^ — k'^ 



In equation f p3|l . the Tan's A-function changes the sign for sufficiently large k and its 
prefactor cannot change the sign. Thus, (k'|u[gf'' |k) is not a non-negative matrix. 



5. Identification of the superfluid transition at the ladder-approximation 
level 

5.1. General description 

We consider the sum of the ladder-type contracted 0-graphs shown in figure HI The 
numerical constant under each term in figure H] is the symmetry number of the 
corresponding ladder-type contracted 0-graph. The algebraic expression of the sum 
of the contracted graphs is given as 

oo 2n 

^i-d = E;r2^ E E ^o(ki)...r/o(k2„) (34) 

n=l ki,...,k2„ o-i,...,<T2„ 

X (ki (Ti , ks 0-2 I Tj^^^ I kg CTg , k4 0-4 )■■■ (k2„_ 1 Cr2n- 1 , k2„Cr2„ I T}^^ I ki (Ti , kg (72 ) • 

We approximate the grand partition function by the ladder-type contracted 0-graph as 

logSy(2;, A/a) ~ -2^Li|(-2;) + Pi^d, (35) 

where Lii(x) := X^^i^"/''^' is the polylogarithm There is no singularity in Li5(— z) 
near the origin and along the positive real axis. To examine a phase transition at the 
ladder- approximation level, we examine a singularity of Vi^d on the real positive axis of 
z. Substituting equation (!26|) into equation we have 

oo 

^lad = - 5Z -^K (Pl,P2)^K (P2,P3) • • -^K (Pn,Pl) , (36) 

n=l K pi,...,p„ 
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^<-'' (2) (8) (24) 

(2" • n) 

Figure 4. Vud as the sum over different ladder-type contracted 0-graphs. The 
corresponding symmetry numbers are shown under the graphs. 

where we define a matrix ^k, related to f/^p and 770 (k), as 

(p, pO := 2V^^^/(^-)(p|<^|pO 

X v/t/o(K/2 + p)r/o(K/2 - p)^r,^{}^/2 + p')r/o(K/2 - p'). (37) 

If IpI , Ip'I ^ A^"*^ In z, (p, p') — 0. Tlien, we can cut off tlie matrices (p, p') 
and consider a sufficiently large square matrix (p, p') for < |p| , |p'| < k'^ 
without any influence on physics by considering of a sufficiently large k'^. The ma- 
trices (p, p') are clearly symmetric (i.e., (p, p') = ^^(p',p)) and primitive 
(i.e., all the elements are non-negative: (p, p') > and its square is positive: 
X]q (Pi q) -^K (Q; P') > 0)- Therefore, the matrices (P? P') satisfy the condi- 
tion for the Perron-Frobenius theorem: 

The Perron-Frobenius theorem 

Let A = (aij) be an n x primitive matrix (i.e., Cjj > and (^'").tj > for some 
natural number m and 1 < i,j < n.). Then, the following statements hold 

(i) There is a positive real number A^^, called the Perron-Frobenius eigenvalue, such 
that A^^ is an eigenvalue of A and the absolute value of any other eigenvalue A 
(possibly, complex) is smaller than A^^, |A| < A^^. Thus, the spectral radius is 
equal to A^^. 

(ii) the Perron-Frobenius eigenvalue A^^ is a simple root of the characteristic 
polynomial of A. Consequently, the eigenspace associated with A^^ is one- 
dimensional. 

(iii) There exists an eigenvector u = (ui,...,u„) of A with eigenvalue A^^, called 
the Perron-Frobenius vector, such that all components of u are positive, i.e., 
Avl = A^^u, ?ii > for 1 < i < n. 

(iv) There are no other positive eigenvectors other than u, i.e. all other eigenvectors 
must have at least one negative or non-real component. 

(v) Let B = ihij) be an n X n primitive matrix and A^^ be the Perron-Frobenius 
eigenvalue of B. If hij > aij for 1 < "i, j < n, then > A^^. 
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Now the matrices A'^ (p, p') can be diagonalized using orthogonal matrices 
f/K(p,p') : 

Uk (p, q) (P> P') Uk (p', q') = Ak (q) 5^,^,. (38) 

p,p' 

Without loss of generality, we assume that the eigenvalue Ak (0) is the Perron- Frobenius 
eigenvalue of the matrix for all K, that is, if q 7^ then Ak (0) > |Ak (q)|- We first 
establish the following lemma: 

Lemma — 

For any q 7^ and K 7^ 0, Ao(0) > |Ak (q)|. 



Proof. Using 770 (k) = ^1 + 2;e ^^2/(2™,)^ ^ have 



^-^KViAm) ^,(K/2 + p)r/o(K/2-p) 

^ j^g/3KV(8m) _^ ^g-/3(p2+K.p)/(2m)j"^ J^g/3K2/(8m) _^ ^g-/3(p2-K.p)/(2m)^ 
(^l + ^g-/3(p2+K.p)/(2m)j~^ ^ ^g-/3(p2-K.p)/(2m) j 



< I 1 + 2e 

-2 



< 1 + ^e-/^P'/(2-) = r/o(p)r/o(-p). (39) 



Hence using equation ( 1371) and inequality (l39l) . we obtain (p, p') < .Aq(p,p'), for 
any p,p'. According to the Perron- Frobenius theorem (v), we have Ak (0) < Ao(0). 
Combined with Ak (0) < Ao (0) and Ak (0) > |Ak (q)|, the lemma is proved. (Q.E.D.) 

Substituting equation (138!) into equation (l36l) . we obtain 

00 ^ 

^^^^ = E-EE(^K(q))^ (40) 

n=l K q 

The radius of convergence of a power series ^1"^=!^^/^ is 1. According to the above 
lemma, if Ao(0) < 1, the power series fHOj) is convergent with the result 

^iad = -EEl°g(l-^K(q)). (41) 

K q 

However, if Ao(0) = 1, then the holomorpy of this function V\a,d is lost, indicating that 
there is a phase transition at Ao(0) = 1. The positive real number Ao(0) changes with 
the parameter of the system such as the temperature, scattering length and number 
density. Therefore, the following theorem is established: 

Theorem 2. — 

Let Ao(0) be the Perron- Frobenius eigenvalue of the primitive matrix 

^^o(p, p') = ^^(p|4?|p') (1 - np(p)) (1 - n^ip')) . (42) 
Then, at Ao(0) = 1, the holomorpy of the function "Piad is lost. 
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According to Theorem 2, solving 

Ao(0) = 1, (43) 

we obtain the critical fugacity = Zc{X/a). We will give an example of this calculation 
in the next section. 

From the Perron- Frobenius theorem (ii) and Theorem 2, it is found that there exists 
a unique fugacity at the phase transition point for any a. Physically this implies that 
there is no fragmentation of Bose- Einstein condensates of dimers and Cooper pairs at the 
ladder-approximation level of Lee- Yang contracted 0-graphs. Besides, we can show [19] 
that the thermodynamic potential, which is calculated at the ladder-approximation level 
of the Lee- Yang contracted 0-graphs and in the first-order approximation in the coupling 
constant of U^'^\ is identical to the results of Nozieres and Schmitt-Rink (NSR) [271128] . 
which include the effect of the Gaussian fluctuations of Cooper pairs. Hence, our result 
indicates that if we only consider the range of Gaussian fluctuations of Cooper pairs, 
there is no fragmentation in the BCS-BEC crossover. 

5.2. Application to Bose-Einstein condensation of strongly hounded dimers (BEC limit) 

We now apply Theorem 2 to the limit of strongly bounded dimers (BEC limit), i.e., 
= w'^ ^ I. We first rewrite equation f l26|) as 

(k'i4f ik) = (k'l^r^ik) + (44) 

Ssc{x) = — - we-'' -^xF (x) + \w\ e"" erfc {\w\) ] , (45) 



(k'|4-''^)|k) = e[w) ■ e-'M^'WM, (46) 

where 6 (x) = {l+x/\x\) /2 and iphO^) = ~~ i+Q^a)'^ • interpret V'b(k) as a 

normalized bound state, because the Fourier transform of iphO^) is the relative wave 
function of the bound state '?/'b(r) = (27ra)~^/^ e~'''/"/r, which is given in equation f l25p . 
Keeping the leading order term, we have 

(k'|4f|k) = (k>g''^|k) + Oiw-') = e-'Mp')Mp) + 0{w-'). (47) 

Then 

^o(p', P) = e"'Vb(p')V'b(p) (1 - nF(pO) (1 - riF(p)) + Oiw~'). (48) 

We obtain 

5^^^(p,q)^b(q) (1 - riF(q)) ^ e"'' ^ [V'b(q) (1 - nF(q))] Vb(p) (1 - n^) ■ (49) 
q q 

Because of the Perron-Frobenius theorem (iii), (iv) and since iphip) (1 — ''^f(p)) > 0, 
iph{p) (1 — ''^f(p)) is the Perron-Frobenius eigenfunction of the matrix Aq. Hence, 
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0.10 0.15 0.20 0.25 0.30 

a/Ac = {'s/^Wc)^^ 

Figure 5. Dependence of = exp{PcfJ'c) on a/A at the zeroth approximation (no 
Fermi statistics, dashed curve) and at the next order (Fermi statistics, soHd curve). 
The critical fugacity is enhanced by the Fermi statistics. 



^2gW)^ X]q[^b(q) (1 — '^f(p))]"' is the Perron-Frobenius eigenvalue Ao(0). Therefore, at 



the phase transition point, 



' e'"^' Yl [^b(q) (1 - n^ip))? = 1 + 0{w^'). 



(50) 



The effect of the Fermi statistics, which is represented by [iphi^) (1 — ''^f(p))] , can 
be evaluated as 



$^[V'b(q)(l-nF(p))]^ = — / dq 



rr' Jo [1 + (qan 



1 + . e^^=^ /(2m) 



n=l 



{1 + nwl) e""""'/^ erfc 



nwt. 



2nw^. 



TT 



1 + 2^/- (Lii(-2c) + Li3(-2c)) -w:-" + 0(10:"). 



(51) 



In the zeroth approximation we do not include any 0{w^ ^) terms. The zeroth-order 



solution for Zr satisfies z"^ e"''^ = 1, and thus we obtain 



l^bl 



(52) 



where -Eb is the binding energy of the dimer. The zeroth approximation thus leads to a 
non-interacting-dimer Bose gas. 

At the next order, we include the terms which contribute to the order of w~^. 
Hence, from equations ( l50l) and f lSTl) we obtain 



2 log z, + wl + 2y^ (hi. {-Zc) + Lia {-zS) ■ + 0{w-^) = 0. (53) 

Here, Lii (— ^c) + Li3 (— 2;^) is negative. Hence, the critical fugacity is enhanced by means 
of the Fermi statistics as shown in figure O 
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6. Summary 

We have applied the Perron-Frobenius theorem to study the superfluid transition of a 
two- component Fermi gas with a zero-range s-wave interaction. We have shown that 
the matrix elements of the two-body cluster functions in momentum representation are 
primitive based on the Bethe-Peierls boundary condition or the Fermi pseudopotential. 
According to the quantum cluster expansion method of Lee and Yang, the grand 
partition function is expressed by the contracted 0-graphs. We have considered 
an infinite series of the ladder-type Lee- Yang contracted 0-graphs and analyzed its 
singularity. It has been found that the singularity is governed by the Perron-Frobenius 
eigenvalue of a certain primitive matrix which is defined in terms of the two-body cluster 
functions and the Fermi distribution functions. As a consequence, there exists a unique 
fugacity at the phase transition point, which indicates that there is no fragmentation 
of Bose-Einstein condensates of dimers and Cooper pairs at the ladder- approximation 
level of Lee- Yang contracted 0-graphs. An application to the EEC of strongly bounded 
dimers has also been made. 
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Appendix A. Derivation of the two-body cluster function for the s-wave 
pseudopotential 

To calculate the two-body cluster function U^'^^ in the s-wave approximation, the general 
formula discussed in Ref. [15] is applied. The formula deals with the case of a central 
potential and the volume V = oo. The relationship between the cases of finite V and 
infinite V is discussed in Ref. [16]. In this section, we show the subscripts V and oo to 
distinguish the cases of finite V and infinite V , respectively. 

The general formula in Ref. [I5] needs the set of energy eigenvalues and 
eigenf unctions, which are calculated by using the Bethe-Peierls boundary condition (l22l) 
or the Fermi pseudopotential (123|1 . Here we use the Fermi pseudopotential ( 123|) . The 
two-body Hamiltonian is 

Hf^ = (V? + V^) + ^6'ir)P. (A.l) 
Here we introduce the center-of-mass and relative coordinates: R = (ri + r2)/2C 
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r = ri — r2, and their absolute values: R :- 
the relative motion is 



|R|, r := |r|. The Schrodinger equation for 



.lv2 + '^S'^{r)^r ) ^(r) = E^(r). 
m m or 



(A.2) 



The solutions of this equation are continuous scattering states ipsi^) as shown in equation 
with energy = k^/m and one bound state ipbii^) as shown in equation f l25l) with 

the binding energy = — l/(ma^). The same results can be obtained by the Bethe- 

Peierls boundary condition. 

Using the general formula for the coordinate representation of the two-body cluster 

function, we obtain 



(r;r'2|t/S,oo|rir2) 



A3 



BP,oo 



where 6 (x) 



C(^)^sc(r') 



x/ /2. The integration over k gives 



sin(/cr) sin(/cr') 



(A.3) 



(A.4) 



(r \u 



(rel) 
BP,oo 



|r) 



Anrr'X 
A 



y2g-m(r+r')V(4/3) 



_|_ ^ g/3/{ma2)g-(r+r')/a g^fp 



r + r' im 



1 I 



(A.5) 



where the complementary error function erfc (x) is defined in equation (130!) . 
The momentum representation can be obtained by using 

1 



d^ridh[dh2dh'2 



(87r3)2_ 

X e-'^'=i('^'--'-'--'^-'-")(r;r'2|t/^'i^,^|rir2), 
where k := (ki — k2)/2, K := ki + k2 and k := |k|. We write 



k'ik'2|t/J,') |kik2) = 6\K - K') (k;k'2|n^^^|kik2) 



'1,'L,(2)| 



(A.6) 



(A.7) 



The function u^"^^ is defined only for the case with k'^ + k2 = ki + k2 and is independent 
of volume. The cluster function defined for a finite volume V is 

r3 



k;k'2|t/^?v|kik2) 



6TT" 



V 



Sk,k' (kik2|n|3^p|kik2). 



We finally obtain 



(k;k'2|ni?) |kik2) = e-^^^/(^-) (kK^lk), 



(rcl) I 



where 



k'|4f |k) 



A s{k') - s{k) 



23/2^^5/2 ^/2 _ j.2 



sik) 



ma 



2ka 



e/3/(™«')erfc 



ma 



(A.8) 

(A.9) 

(A.IO) 
(A.ll) 
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which gives equation 

The two-body cluster function [/gp ^ for the s-wave pseudopotential has been given 
in Ref. [17] for positive a, but the resuh of this Appendix holds for arbitrary a. 

Appendix B. Proof of the strictly increasing of the function 

In this Appendix, we prove that ^s(x,w) > for < x < oo and — oo < w, where 
^s{x,w) is the derivative of the equation (!28|) which is written as 



d , . 2 X 



-S X, w 



1 



i-w) e"^ (1 + + w^) - 1 - -F (x 



2xF (x) (1 + x^ + w^) - ; \ e"''erfc (-w) . (B.l) 



x^ + 

We will discuss separately the following four cases: (i) — oo < w < —1, (ii) —l<w< 
-|, (iii) -| < w < 0, (iv) 0<w. 

We first organize inequalities which are related to the Dawson's integral and the 
complementary error function. The lower and upper bounds for the Dawson's integral 
are derived, by considering the Taylor series of the integrand of the Dawson's integral 
and estimating each term, as 



1 1 e~^^ /x^ 1 \ 

, , 1 1 3 3e-^' / 5x2 5^6 \ 

F (x) < e"^' (^x+'^x^] (0 < X < 1) 

F (x) > e'""^ ^x + ^x^ 

The upper bounds for the complementary error function are derived as follows: 

1 e"^"" 

erfc(x) < —p= (0 < x), 

Vvr X 

erfc(x) < ^e"""^ + ^e"^"^^ < x 

2 1 _ 2 
erfc(x) < —= e""" (0 < x), 

VTTX + 1 



B.2 
B.3 
B.4 
B.5 

B.6 

B.r 

B.8 

2 2 

erfc (x) < 1 ^x H ^x^ (0 < x) . (B.9 

(i) —00 < w < — 1 

The inequalities (EJ), f lB^l) and fiRBll lead to 

^ 1 X ~ 

— s(x,m;) > — ^^i(x,w), (B.IO) 

ox yJH (X^ + W^Y 
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where 
S'i(x, w) 



w 

-2 W - 2y/n(l + + w^)w H \ . 

8 4 4 



For — oo < w < — 1, we notice that 



^S,{x,-l) = -6^+8 
ow 



12^/^w - 8 + — > 12^/^ - 8 > 0, 



20F(1 + x^) < + 8 < 0, 



(B.ll) 

(B.12) 
(B.13) 



S,{x, -1) = -2 - ^a;2 - 4 + 2v/^(2 + a;^) + ^ + ^ 
8 4 4 

= 4v/^-6+ (^2v^+^-^^x2>0. 
Thus, ^[(x,!^) > 0. Therefore, if — oo < w < —1, then ^s{x,w) > 0. 
(ii) -1 < w < -i 

The inequahties (ES]), f lB^l) and fiRTll lead to 



(B.14) 



d , , 1 
— > 
ox 



71 xiX'^ + W 



2\2 



Sn{x,w), 



(B.15) 



where 



S'ii(x, w) = + wh2{w)x + w ( hs{w) 



+2x^ 



(l + y/Tlwhi{w)) 



1 



1 — 

1 + If 



2-/Tihi{w) ) a; 



1 — 



-X' 



= e"' erfc(— w) 

h2{w) 

1 

4 



4 4w ^ 
1 — u; 3 



1 — w 



2y/Ti - 2^w^ -4w;. 



(B.16) 

(B.17) 
(B.18) 

(B.19) 



77^ + ^ > 

\\~wY 3 



We will evaluate each term in equation (lB.16p . We notice that hl^iw 
and /i2 (-|) = 2 - 2^ < 0. Thus, h2{w) < for < -\. The inequality (IRT]) lead 
to h^{w) < I + le-s*"" for w < -\. Thus 

/i3(w) 



^ 2v/7r/ii(uO < -4w- 2y/TTw'^ + A/vre""^ = /i4(w). 



1-w -y-'- 3 

We notice that {-^ < 0, h'^ (-i) > and h'l{w) = -Ay/l^-^e-"^ (l - ^) < 



(B.20) 




for —l<w< — ^. Thus, h^iyo) < holds. Moreover, the following inequality holds: 



1 + a/vtw hiiw) > 



1 + w 
1 — w 



{w < 0). 



(B.21) 



Hence, S\\{x^w) > 0. Therefore, if — 1 < w < — i, then -^s{x,w) > 0. 
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(iii) -\<w <Q 

For — i < < 0, we put —w = ax {a > 0) and discuss separately the following three 
cases: (iii-a) | < x, (iii-b) < x < | and 1 < a, (iii-c) < x < | and < a < 1. 



(iii-a) We have < a < 1. The inequalities (ESI), (El) and flRSl) lead to 



d 1 e"^ 1 

dx ' ~ ^/tt x^ {1 + a'^y 1 + ax 



Sm{x,a), 



where 



S]]](x, a) = a 



1 



x' + 



9 4a2 



4 



2y^a^{l + a^ 



2(2 — A/7r)Q;(l + a'^)x^ — 2(2 - v^)aV - 2(2 - ^/^) 



X 



ax 



If I < X and < a < 1, then 



dx' 



■Siii{x,a) > 



135 , 

— - 22 2 - a 
lb 



X > 0, 



(B.22) 



(B.23) 



(B.24) 



d . (I 



3v^ 



-a 



53\ 3 fll^ 7 



TT a + , 

487 V 8 



17 r- 108n 27 ^ 

— VTT a H > 0, 

8 ^ 256 / 64 



(B.25) 



32 



-a 



+ 



/50r _ 319\ 
V^6~ "5127 



v^_ll-^ 
16 96 
9 



/5v^_ 11 
V 32 48 



a 



a 



256 



> 0. 



(B.26) 



Here, inequalities (IB.25P and f lB.26|) . in which not all coefficients are positive, follow 
from evaluating the cubic functions for < a < 1. Hence, S'iii(x, a) > 0. Therefore, if 



x>\ and — I < w < 0, then ^s(x, w) > 0. 



(iii-b) For -i < < and < x < i, the inequalities (IB.4p . flB.5p . flB.9P and 

^a^x^ < 1 + a^x^ + a'^x'^ (for < ax = —w < 1) lead to 



^s{x w) > ^ 

dx ' ~ (1 + a^)2x^ 

2x3 



a H — ^(1 — a 

TT 



ax 



30F L 



(3 + 2a2) e-" + 2a^ 



- a^x^ + 



2x^ 



(1 + a^)e 



(1 + a')e-^' + 2a6 - aV }. (B.27) 



a 



Now, we consider the case of 1 < a. Using flB.27p . e ^ — 1 > — x^ and x < l/(2a), we 
obtain 



d . . 4 X 
ox 3v7r 1 + a^ 

Since 

3 + 2a2 
1 + a2 



3 + 2a^ 
l + a2 



-l/(4a2) _ 



3v^ 



3v^\ 



a 



4 / l + a2_ 



5 
-( 
2 



3v^ 



(B.28) 



(B.29) 
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and 2 > 3-\/7r/4, we obtain ^s{x,w) > 0. 



(iii-c) The inequality ( IB. 271) obviously holds. Using ( 1B.27I) and e~^' > 1 — > |, we 
obtain 

X 



dx'^'"''"^- v/^(l + «2)2 



[W + 1 - ^/Tia{a^ + 1 



If < tt < 1, then + 1 - ^a{a^ + 1) > 0.7 > 0. Thus, j^s{x, w) > 0. 
(iv) 0<w 



d_ 

dx ' 



(x, w] 


as £s( 




d ^ 
dx 


[x,w) = 


2w X 


2^ 


d 

OX 


[x,w) = 


2.^ r 









dx 

2,^2 



w +x 



- (1+W'' + X^ 



-e'^'erf (w) - -^F(x) 

71 



X 



[{2x^ + l)F{x) -x] 



Here erf {x) := j^dte *^ = 1 — erfc(x) is the error function. 
The following inequalities hold: 

> 1 + ^ + 
e"'"erf iw) > 



, erf (w)>l> 
2w ^ ^ X 

„2 



(2x^ + l)F(x) -X > (x>0). 



(B.30) 



(B.31) 



(B.32) 



(B.33) 
(B.34) 
(B.35) 



Thus, for < w, we obtain -^si{x,w) > and 4-S2{x,w) > 0. Then -^s{x,w) > 0. 
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